Strain mediated suppression of the metal-insulator transition in EuNi0 3 thin films 
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Ultrathin epitaxial films of EuNiOs were grown on a series of substrates traversing highly compressive (- 
2.4%) to highly tensile (2.5%) lattice mismatch. X-ray diffraction measurements showed the expected c-lattice 
parameter shift for compressive strain, but no detectable shift for tensilely strained substrates, while reciprocal 
space mapping confirmed the tensile strained film maintained epitaxial coherence. Transport measurements 
showed a successively (from tensile to compressive) lower resistance and a complete suppression of the metal- 
insulator transition at highly compressive lattice mismatch. Corroborating these findings, X-ray absorption 
spectroscopy measurements revealed a strong multiplet splitting in the tensile samples that progressively weak- 
ens with increasing compressive strain that, combined with cluster calculations, showed enhanced covalence 
between Ni-d and O-p orbitals leads to the metallic state. 



I. INTRODUCTION 



Complex transition metal oxides are at the forefront of cur- 
rent research towards the understanding of the fundamental 
physics underlying several astonishing physical properties in- 
cluding high temperature superconductivity, colossal magne- 
toresistance, multiferroicity, and the thermally induced metal- 
insulator transition (MIT)[ - ]. In particular, the tempera- 
ture driven MIT in correlated oxides has garnered remark- 
able investigation over last several decades. Understanding 
of the MIT and its control by external means such as pres- 
sure, magnetic field, confinement, and chemical doping is not 
only interesting from the physics point of view[7-10], but 
they also show great potential as candidates for future elec- 
tronic devices [ ]. Recent advances in material synthesis by 
utilizing the lattice mismatch between the substrate and the 
as-grown films have opened a new dimension in controlling 
and engineering the materials properties. This epitaxial strain 
leads to a restructuring of the grown material and an entirely 
coherently strained crystal can be obtained in the ultra thin 
limit. The effect of this lattice modulation on the materi- 
als properties can be quite dramatic[ 12-25] and is of partic- 
ular interest as technology continues to accelerate towards the 
need for ultrathin film materials for electronic applications, 
and, thus, the need for a fundamental understanding of these 
effects is paramount. 

The rare-earth (RE) nickelates ReNi0 3 (Re = Pr, Nd, Eu, 
...) in their bulk form, with Ni having formal +3 oxidation 
state (t^Le 1 ), all display a metal-insulator transition, except 
Re = La, while the nature of the transition and its tempera- 
ture (Tmit) depends strongly on the choice of Re cation as 
shown in Fig. 1 [26-28]. The very distorted ReNiC>3 with Re= 
Lu, Y, Eu, Sm, etc. first show a second order MIT at higher 
temperature while the spins remain in their paramagnetic state 
across the transition. These compounds further undergo an- 
other second order transition to a E' - type antiferromagnetic, 
insulating state at lower T (Tjy) along with possibly forming a 
charge ordered state, which is still controversial and not fully 
understood[29-33]. On the other hand, the members with in- 
termediate distortion {e.g. Nd, Pr, Smo.sNdo.s) exhibit a first 
order phase transition, accompanied by a small volume expan- 



sion, directly from the paramagnetic metallic state to the E' - 
type antiferromagnetic insulating state, bypassing entirely the 
paramagnetic insulatoring region[27, 34-36]. Several tech- 
nological applications have already been proposed and de- 
vices fabricated utilizing thin films of these materials [37-40]. 
From a technological point of view, control of the MIT is of 
great importance, however using different RE ions depending 
on application is impractical, as synthesis conditions and ther- 
mal stability can vary among them[12, 41^-4]. Alternatively, 
it would be more practical, as well as more scientifically ap- 
proachable, to use this obligatory strain in thin films as an 
avenue to engineer the physical properties within a particu- 
lar RE nickelate as it can modulate the MIT and AFM transi- 
tion significantly as observed for NdNiC>3 [12, 13]. However, 
several key questions still need to be answered. Key among 
them, what effect will strain (tensile and compressive) have 
upon the phase diagram in Fig. 1 for the materials showing 
second order transitions? It has been shown, through powder 
studies, that external pressure suppresses the MIT while ac- 
tually slightly raising the AFM transition, eventually leading 
to a surprising AFM, metallic state[45-47]. It is also possible 
that epitaxial strain can be used to take one of the more dis- 
torted nickelates, showing two second order transitions, and 
induce a first order phase transition via strain mimicking the 
results of chemical doping[48]. 



In this letter, we characterize several ultra thin epitaxi- 
ally strained EuNiC>3 (ENO) films grown on various sub- 
strates spanning both compressive and tensile strain using x- 
ray diffraction (XRD), reciprocal space mapping (RSM), elec- 
tronic transport, and x-ray absorption spectroscopy (XAS) to 
elucidate the effects of epitaxial strain on the structural, elec- 
tronic, and magnetic properties. The results of these methods 
display the capacity of epitaxial strain to modulate the physi- 
cal properties of this highly distorted nickelate system. In ad- 
dition, cluster calculations on a NiOg octahedra allowed the 
extraction of the charge-transfer (CT) energy on the various 
films, shedding light on the source of the observed modula- 
tion in material's properties. 



II. EXPERIMENTAL 

ENO films were grown on a variety of substrates incor- 
porating lattice mismatches ranging from 2.5% to -2.4% us- 
ing pulsed laser interval deposition, the details of which can 
be found elsewhere [44]. The substrates used for growth 
presented here are as follows: YAIO3 (YAO; -2.4% lat- 
tice mismatch), SrLaA10 4 (SLAO; -1.3%), LaA10 3 (LAO; 
-0.3%), NdGa0 3 (NGO; 1.5%), and LaGa0 3 (LGO; 2.5%). 
XRD measurements, besides ENO on YAO which is de- 
scribed elsewhere[44], were taken with a Panalytical X'Pert 
Pro MRD (Panalytical, Almelo, The Netherlands), equipped 
with a parabolic mirror and triple bounce / axis monochro- 
mator on the incident and diffracted beams, around the (002) 
(psuedocubic notation) truncation rod of the substrate . The 
same instrument was used to measure a RSM around the (- 
103) truncation rod. Transport properties were measured with 
a Quantum Design Physical Property Measurement System 
(PPMS) using a four point probe with the Van Der Pauw ge- 
ometry. XAS measurements were taken at the 4-ID-C beam 
line of the Advanced Photon source at Argonne National Labs 
in total electron yield (TEY) mode on the Ni L2.3-edges at 
250° K. 

Calculations were done for a NiOg cluster with octahe- 
dral coordination using the methods described here[ ]. The 
Hamiltonian includes the on-site Coulomb interaction be- 
tween the 3d electrons and between the 3d electrons and the 2p 
core hole. The parameters are calculated within the Hartree- 
Fock limit and scaled down to 80% to account for intra-atomic 
screening effects. The monopole parts were F dd = 6 and ¥ Q d = 
7 eV. The spin-orbit coupling was included for the 3d and 2p 
electrons. The hybridization with the ligands was taken into 
account by including configurations up to double ligand hole. 
The hybridization parameter is V = 2:25; -1:03 eV for the e g 
and t 2ff orbitals, respectively. The cubic crystal field of 10 Dq 
is 1.5 eV. 



III. RESULTS 

Fig. 2(a) shows 26 - oj scans around the (002) truncation 
rod for the 15 unit cell (uc) ENO films grown on different 
substrates ((006) rod for SLAO, due to its tetragonal struc- 
ture). All films show a broadened film peak, due to the 15 uc 
thickness of the films, and a sharp substrate peak which was 
used to align each data set. The film peaks for the highly com- 
pressive films (YAO, SLAO) show a noticeable shift from the 
ENO lattice constant (represented by a dashed line) towards 
smaller 26, while for LAO no shift is observable. Films un- 
der tensile strain (NGO, LGO) show a well resolved film peak 
with no shift from the bulk value. Fig. 2(b) shows a RSM 
around the (-103) Bragg peak for a 35 uc ENO on NGO film. 
The weak film peak shares the same value for H (reciprocal 
lattice units) and a larger value for L. 

DC transport was measured on the 15 uc films, shown in 
Fig. 3(a). The data was recorded during both cooling and 
heating cycles from 380 K to 2K, though no hysteresis was 
found and only warming curves are shown. For the tensile 



cases the resistivity follows the expected bulk-like insulating 
behavior below 380K. For the small compressive strain case, 
LAO, the resistivity at lower temperatures separates from bulk 
behavior ~ 250K and begins increasing at a lesser rate. For 
the compressively strained SLAO the sample shows metal- 
lic behavior at high temperatures with a MIT occurring at 
335K. The film on YAO is metallic down to 2K. Fig. 3(b) 
vs T for all films besides ENO on YAO, with each 



shows 



dlnp 



showing a characteristic kink indicative of an AFM transition 
around 200K. Fig. 3(c) shows the extracted T values for the 
kinks, denoted T*. 

XAS measurements were taken on the Ni L3.2 - edges, plot- 
ted in Fig. 4(a). All films show a strong white line at ~ 855 
eV (~ 872 eV). Additionally, a shoulder around 853 eV (871 
eV) is apparent for all films, being much larger for the ten- 
sile case and gradually decreasing as compressive strain is in- 
creased. A small peak around ^851 eV corresponding to the 
La M4 edge appears for the substrates containing La. The size 
of the energy speration between the L3 and multiplet peak for 
each strain is shown in Fig. 4(b) in the left axis. The splitting 
decreases from the tensilely strained films, ^1.8 eV, with in- 
creasing compressive strain to ~ 1.2 eV for YAO. Fig. 4(b) 
also shows the calculated CT energy (right axis) which fol- 
lows a very similar trend. 



IV. DISCUSSION 

Epitaxial strain is possible with these materials due to the 
extraordinary ability of the perovskite lattice to accommodate 
the strain through lattice and symmetry effects; it is these ef- 
fects that ultimately lead to the observed modulation of the 
physical properties. Thus XRD is a very important technique 
to deduce how the lattice incorporates the strain and to check 
for secondary phase formation. 26 values, corresponding to 
the film peaks in Fig. 2(a), were used to calculate the out- 
of-plane lattice constants. This yielded the following c-axis 
lattice parameters: 3.878 A (YAO), 3.84 A (SLAO), 3.80 A 
(NGO), and 3.81 A (LGO), while for LAO strong overlap 
between the strong substrate peak and film peak due to the 
small strain value, -0.5%, prevent a reliable c-axis lattice con- 
stant from being extracted and the weakness of the film peak 
is probably due to the severe twinning present on the LAO 
substrates. While for the samples with in-plane compressive 
strain the shift of the out-of-plane lattice constants from the 
bulk value of 3.80 A is expected and consistent, the tensilely 
strained substrates show no significant shift (0.26% for ENO 
on LGO, much lower then the 2.5% biaxial in-plane strain). 

In order to ascertain whether or not this lack of c-axis lattice 
modulation was due to strain relaxation, unlikely in such thin 
films, a RSM was taken for a 35uc ENO film grown on NGO, 
Fig. 2(b). The thicker sample was needed in order to obtain a 
strong enough (-103) film peak with this in-house MRD sys- 
tem, and, being thicker than 15 uc, was more likely to show 
relaxation. The RSM around the (-103) reflection of the sub- 
strate shows a strong substrate peak with a broadened film 
peak at higher L (shown in reciprocal lattice units (r.l.u.)). The 
H value of this film peak matches well with that of the sub- 



strate (demonstrated by the dotted line) confirming the film is 
coherent to the substrate, while the center L value of ~ 3.053 
r.l.u. gives an out-of-plane lattice constant of 3.793 A, in good 
agreement with the rocking curve measurement. This, along 
with the bulk-like c-axis lattice constant, implies the strain is 
compensated for by some other means, possibly as was found 
for LaNiOa films under tensile strain[50]. In that case the oc- 
currence of "breathing" modes were found that compensated 
for the a-b lattice mismatch without affecting the c-axis lat- 
tice constant. X-ray linear dichroism (XLD) measurements 
and density functional theory calculations will be done in the 
future to confirm this. 

Tracking the evolution of the MIT with lattice modulation 
in these epitaxially strained films revealed a very significant 
effect. As compressive strain is increased Tmit is suppressed 
until entirely disappearing for the case of ENO on YAO (seen 
only for Re = La in bulk), highlighting the ability of utiliz- 
ing lattice mismatch to engineer the behavior of physical phe- 
nomena. This sample was found to have linear T dependence 
above 200K, which is common to the ReNiC>3 compounds in 
their metallic phase excluding Re = La, followed by T 3 ' 2 de- 
pendence below. The T 3 / 2 dependence at low T was shown to 
be caused by enhanced scattering due to bond-length fluctu- 
ations, indicating strong electronic correlations and enhanced 
vibronic fluctuations of the charge ordered state, likely due to 
quantum confinement in these ultra thin films[5 1, 52]. In the 
case of SLAO, with intermediate compressive strain, linear T- 
dependent metallic behavior is followed by a MIT shifted to 
335K, putting it enticingly close to room temperature which 
makes it a potential candidate for electronics applications. In- 
terestingly, the tensile strained samples have nearly identical 
resistivities for all temperatures observed despite their con- 
siderably disparate strain values, possibly associated with the 
lack of c-axis deformation as observed with XRD. Unfortu- 
nately, the high temperature of the bulk MIT (480 K) prevents 
us from investigating the change in the Tmit for LAO, NGO, 
and LGO films. With the lack of hysteresis, characteristic of 
first order phase transitions, our results suggest epitaxial strain 
does not induce a first order transition in this material, as was 
seen with chemical doping[ ]. Instead they show compres- 
sive strain works to lower the transition temperatures as was 
seen with external pressure[35, 46, 47]. This is also strongly 
reminiscent of work on ultra thin films of NdNK>3, where it 
was proposed that a closing of the CT gap was responsible for 
the quenching of the MIT with strain, analogous to changing 
the Re ion in bulk[ 1 2, 27]. To further investigate the possibil- 
ity of epitaxial strain mimicking external pressure the AFM 
transition temperature, T^r, is needed.. 

In an attempt to locate the AFM transition temperature, 



was investigated for all insulating samples, shown in 



dlnp 
d(J) 

Fig. 3(b). This analysis helps reveal the temperature where 
spin ordering occurs, and appears as a small kink in the 
T dependence [ ]. Our films all exhibit a broadened kink 
around T*, shown in Fig. 3(c). The change in temperature 
of this feature was very small, approximately 12K lower for 
the tensile strained samples, similar to the changes reported 
for external pressure [ 5]. For the tensilly strained samples a 
lower value of 190K was found, possibly due to the thin film 



nature of the films[53]. For LAO, with very small compres- 
sive strain, the value is shifted upward to 203 K. For SLAO the 
value is shifted to 207 K, strongly suggesting the Neel tem- 
perature for these films is shifted higher with higher strain, 
analogous to the effect of external pressure found by Zhou 
et al, however a different mechanism must be responsible, as 
A is not constant, discussed below. Resonant x-ray scatter- 
ing measurements are needed in order to further confirm these 
changes in T^v and investigate any possible changes in mag- 
netic structure, which cannot be investigated via transport. 

XAS is an excellent experimental tool to study electronic 
structure which has been extensively utilized in the study of 
ReNiOs perovskites[ 1 2, 54, 55]. Though TEY mode is a sur- 
face sensitive tool, the smaller thickness (~ 6 nm) of these 
films compared to probing depth (~ 12 nm) of TEY mode al- 
lowed us to explore the electronic structure of the entire film. 
The strong white line at 855 eV (872eV) corresponds to the 
L3 (L2) edge transition from the d 7 : tfLeJ ground state to the 
ctjLe 2 , (c denotes a core hole) excited state. Another lower 
energy peak, ~ 853 eV, which corresponds to the same elec- 
tronic transition, is due to a multiplet interaction between the 
core hole and the excited electron when strong localization is 
present, becoming prominent for the highly tensile strained 
samples. This systematic change in the multiplet / L3 rela- 
tive position and intensity holds valuable information about 
the hybridization of the d 7 and d 8 L states. The quantitative 
value of the observed splitting, which is found to be related 
with the CT energy (A) has been plotted in Fig. 4(b)[ ]. 
As can be seen, this splitting begins decreasing as the films 
are further compressed in the a-b plane, indicating a change 
in the covalence between Ni and O is likely causing the ob- 
served transport properties. To confirm this explanation cluster 
calculations were carried out using A as a a control parame- 
ter, Fig. 4(b) (right side). As shown changing A reproduces 
the observed splitting, with a corresponding value of ~ 2.5 eV 
(for the tensile strained samples) and being reduced down to 
~ 0.65 eV in the case of YAO. 

The large reduction in A (approximately l/4th of the sat- 
urated value for ENO on YAO) implies the covalence is 
strongly increased between Ni-d and O-p orbitals. This can 
be understood in terms of changes in the Ni-O-Ni bond; as 
the lattice is compressed in plane an increase in the overlap 
of the O: p XtV and Ni: d I ,2_ y 2 orbitals occurs. The increasing 
covalence is strongly reminiscent of what was seen for A-site 
cation exchange (Fig. 1)[ ]. For the tensilely strained sam- 
ples, using similar analysis, the overlap would decrease, lead- 
ing to a reduction in covalence, however the splitting appears 
to remain unchanged with increasing tensile strain. These re- 
sults indicate that as the strain is changed from compressive to 
tensile the d 8 L state is effectively decoupled from the d 7 state. 
Park et al have recently proposed, using DFT + DMFT, a site- 
selective Mott state, where the insulating gap is determined 
by the singlet formation energy between an O-p hole and Ni-d 
electron[ ], our results are compatible with their predictions, 
showing the decoupling of the d 7 and the d 8 L states leads to 
an insulating state. Further corroborating our results, Wang et 
al recently showed that the insulating regime is largely con- 
trolled by the number of d electrons, not the intraelectronic 



repulsion, which is strongly dependent on the charge transfer 
from oxygen atoms [ ]. To further investigate this effect and 
confirm our conclusions, resonant x-ray scattering measure- 
ments will be performed to investigate the evolution of the 
bond lengths and angles with strain. 



V. CONCLUSION 

To summarize, epitaxial ultrathin films of ENO grown on 
a variety of substrates spanning both compressive and tensile 
strain were investigated with XRD, electronic transport, XAS, 
and cluster calculations. The sample's properties were found 
to be highly tunable with strain, shifting the MIT tantalizingly 
close to room temperature and only effecting Tjv to a small 
extent. These results indicate the transitions remains second- 



order, showing a key difference between epitaxial strain and 
A-site exchange or doping and implying the application of 
compressive strain mimics external pressure. Through XAS 
and cluster calculations it was determined that compressive 
strain enhances the covalence of the Ni-d and O-p orbitals, 
leading to an entirely metallic state. These results showcase 
ENO's potential utility for future device applications, as the 
MIT can be shifted from the high bulk value of 480K to near 
room temperature or to being entirely quenched without the 
complication of using different materials or chemical doping. 
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FIG. 1: (Color online) Phase diagram for the nickelates of different 
A-sites, data taken from Lengsdorg et al[ ] 
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FIG. 2: (Color online) (a) XRD data for 15uc ENO samples on vari- 
ous substrates. The shifting from the bulk lattice value (indicated by 
the dashed line) is apparent for the compressively strained samples, 
(b) RSM for a 35uc ENO film grown on NGO showing the film is 
coherently strained. 
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FIG. 3: (Color online) (a) Transport data for 15uc ENO samples on 
various substrates. The dashed line indicates the temperature where 
the XAS data was taken (Fig. 4a). (b) dln(p)/d(l/T) data for the 
films, (c) T* for various strains. 
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FIG. 4: (Color Online) (a) XAS data for 15uc ENO samples on var- 
ious substrates showing the change in the pre edge peak with strain, 
(b) The experimentally obtained peak splitting along with the theo- 
retically obtained corresponding CT energy. 



